THz lasing in a polariton system: Quantum theory 
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We study the laser regime of terahertz (THz) emission from a semiconductor microcavity in the strong cou- 
pling regime, where optical transitions between upper and lower exciton-polariton modes are allowed due to the 
mixing of the upper mode with one of the dark exciton states. Using a system of master-Boltzmann equations 
describing both polariton modes and the THz mode, we calculate the first and second order coherences and the 
spectral shape of THz emission. This analysis shows that THz lasing in microcavities is possible provided that 
the system is embedded in a good THz cavity, that the optical (polariton) lasing condition is fulfilled and if the 
depletion of the upper polariton mode due to acoustic phonon assisted relaxation processes is reduced. This 
latter condition is likely to be realised in pillar microcavities, which seem to be the most suitable candidates for 
realisation of a THz laser. 



Realisation of efficient terahertz (THz) radiation sources 
and detectors is one of the important objectives of modern 
applied optics (TJ |2). THz radiation based methods have a 
lot of potential applications in biology, medicine, security and 
non-destructive in-depth imaging. Also, wireless data transfer 
utilising THz radiation could provide higher transfer rates for 
in-door short distance or high altitude communications. 

To allow such applications of THz radiation, the creation 
of cheap, reliable, scalable and portable emitters is extremely 
important. None of the existing ones so far satisfies all the 
application requirements. For example, the emitters based on 
nonlinear-optical frequency down-conversion, gas THz laser, 
vacuum tube and systems based on short-pulse spectroscopy 
are bulky, expensive, and power consuming. Various semicon- 
ductor devices based on intersubband optical transitions 
are compact but have a limited wavelength adjustment range, 
have low quantum efficiency and require cryogenic cooling. 
Among the factors which limit the efficiency of semiconduc- 
tor THz sources is the short lifetime of the electronic states in- 
volved (typically, fractions of a nanosecond) compared to the 
time for spontaneous emission of a THz photon (typically mil- 
liseconds). The methods of reducing this mismatch include 
the use of the Purcell effect in THz cavities or the cascade 
effect in quantum cascade lasers (4j-[6l (QCL). Nevertheless, 
till now the QCL in the spectral region about ITHz remains 
costly and short-lived and still show a quantum efficiency of 
less than 1%. Moreover, so far there are no commercially 
available reliable compact and cheap detectors of THz radia- 
tion which are in great need for information communication 
technologies. 

Recent studies of strong coupling intersubband microcav- 
ities have shown the possibility of stimulated scattering of 
intersubband polaritons [7]. Very recently, it was proposed 
to generate THz radiation in semiconductor microcavities in 
the regime of exciton-polariton lasing [8]. The quantum ef- 
ficiency of this THz source is governed by population of the 
final polariton state, which may be tuned over a large range by 
means of the optical pumping. In the strong coupling regime 
in a microcavity 0, the dispersion of exciton-polaritons is 
described by two bands both having minima at zero in-plane 
wave vector k. At k = 0, the energy splitting between the 
two branches is of the order of several meV, which makes this 
system attractive for THz applications (ITHz corresponds to 




FIG. 1. Scheme for the THz laser in a microcavity embedded in 
a THz cavity. Upper and lower polariton branches (dispersion rela- 
tion) are depicted in blue and pink, respectively, with the mechanisms 
considered for all the modes (a, b, c)\ decay (y, y c ), non-resonant 
continuous pump (electronic injection) on the upper branch (P) and 
upper to lower polariton relaxation via exciton reservoir (y a b)- THz 
photons are emitted (absorbed) during the upper to lower (lower to 
upper) polariton conversion at rate G. 



4meV). Stimulated scattering of exciton polaritons into the 
lowest energy state leads to so-called polariton lasing, re- 
cently observed in GaAs ifTOl and GaN ifTTIl based microcav- 
ities. If the scattering from the upper to the lower polariton 
branch were accompanied by the emission a photon, polari- 
ton lasers would emit THz radiation, and this emission would 
be stimulated by the population of the lowest energy polariton 
state. However, this process is a priori forbidden since an opti- 
cal dipole operator cannot directly couple the polariton states 
formed by the same exciton state. This obstacle can be re- 
moved if one of the polariton states of interest is mixed with an 
exciton state of a different parity, say, the elhh2 exciton state 
formed by an electron at the lowest energy level in a quantum 
well (QW) and a heavy hole at the second energy level in the 
QW 0. This state is typically a few meV above the exciton 
ground state, elhhl. Nevertheless, by an appropriate choice 
of the QW width and exciton-photon detuning in the micro- 
cavity, the state can be brought into resonance with the lowest 
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energy upper polariton state. Being resonant, the two states 
can be easily hybridised by any weak perturbation, such as, 
e.g., a built-in or applied electric field. The optical transition 
between such a hybridised state and the lowest elhhl exciton 
polariton state is allowed. 

In this Letter, we describe theoretically the quantum prop- 
erties of THz light generated by polariton lasers according to 
the above scheme. We consider the model system shown in 
Fig. [T] It consists of a microcavity with embedded quantum 
wells operating in the strong coupling regime, and placed in- 
side a THz cavity. Together with the wave-guiding effect of 
the microcavity structure, adding the lateral THz cavity would 
achieve an effective 3D confinement of the THz mode, giving 
rise to enhancement of spontaneous emission rate through the 
Purcell effect U2HII). We will show that in the polariton las- 
ing regime, this device also operates as a THz laser. Moreover, 
polariton lasing may be triggered by an external THz radia- 
tion which would stimulate scattering of polaritons between 
the upper and lower modes in the microcavity. Therefore, the 
proposed device could also operate as a detector of THz radi- 
ation. 

The two polariton branches at k = are modelled by 
two Bose fields with annihilation operators a (upper) and b 
(lower). Both modes (being at resonance) lose particles at 
the same rate 7. This rate provides the units of the problem 
(for standard microcavities, 7^ O.lps -1 ). In Fig.[l] the sys- 
tem is sketched together with the relevant parameters. Only 
one of the modes, a, receives incoherently polaritons at a rate 
P due to electronic injection. Upper polaritons can convert 
into lower polaritons by emitting a THz photon (with oper- 
ator c) at a rate G. Energy is conserved in this process as 
E c = E a — Eh. The opposite process happens with the same 
probability G since the THz photons, once emitted, remain in 
the cavity long enough to make the transformation reversible. 
The THz cavity is not perfect and does, however, lose particles 
at a rate y c . There is a last relaxation process that must be in- 
cluded: the conversion of upper polaritons into lower ones via 
the exciton reservoir at rate 7^, through irreversible phonon 
emission. This process is important in planar microcavities 
and is unfavorable for THz generation, as will be shown be- 
low. Fortunately, the acoustic phonon induced depletion of the 
upper polariton mode may be strongly suppressed in micropil- 
lars |[T0| . All these incoherent processes can be described as 
Lindblad terms (jSf) in the master equation of the total density 
matrix of the system, p : 
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where Jzf^p = 20pO^ - O^Op - pO^O. All the irreversible 
processes are gathered in Eq. ([T^), that is, decay, non-resonant 
pump and relaxation processes. Eq. ([TJ)) represents the non- 
linear processes of upper to lower polariton conversion and 
THz photon emission (and vice versa) similarly to the evapo- 
rative cooling processes in an atom laser fT5lfT7l . 

The Hilbert space of the three modes is spanned by the 
number states { | n, m, r) with n, m, r G N} the number of upper 




FIG. 2. Relevant quantities in the steady state as a function of pump, 
for the THz field (x = c in blue, circles), the upper polaritons (x = a, 
in purple, squares) and the lower polaritons (x = b in brown, rhom- 
bus), (a) Average populations n x . (b) Second order coherence func- 

(2) (2) (2) (2) 

tions gx ' . (c) Cross correlation functions g bc J , g K ac J and g K ab J . (d) 
Linewidth of the THz spectrum of emission. Parameters are: G = 7, 
Yab = and y c = 0.1y. 



polaritons, lower polaritons and THz photons, respectively. 
The occupation of these states described by the probability 
function ^[ njm?r ] = (n,m,r \ p \ n^m, r) follows the set of cou- 
pled quantum Boltzmann master equations filfl : 

dt ^[n : m,r] = - { j{n + m) + y c r + P(n + 1 ) } &^ mA 
+ r{ {n + 1) + (m + 1) ^, m+ i, r] } 

+ Yc (r + 1 ) &[nm+l] + Pn &[n- 1 ,m,r] 
+ Jab { + l)^^Vl,m-l,r] ~n{m+ 1) } 



r(n-\- l)mr|. 



[/2+l,m— l,r— 1] 



■9>\ 



n,m,r] J* 



+ Gn(m+ l)(r+ l){^_i, m+ i, r+ i] - ^ [w] } • (2) 

In what follows, we solve these equations numerically 
in the steady state of the system. Then, any average 
single-time quantity can be computed as (O) = Tr(pO) = 
Hn,m,r ^[rc,m,r] m ? r I O | ^,m,r). The relevant quantities 
characterizing each mode (x = a,b,c), that we use to study 
lasing emission and correlations are: the average occupation 
numbers n x = (xtx); the second order temporal correlation 

function at zero delay, g x 2 ^ = (x^x^xx) /n*, that is 2 for a ther- 
mal and 1 for a Poissonian distribution; the cross correlation 
function between two modes x and y, g$ = (x^xy^y) / (n x n y ). 
These quantities are plotted in Fig. [2ja,b,c) as a function of 
the pumping rate P. Let us explain the choice of parameters 
of this figure, that brings the system in the THz lasing regime. 

The first condition for THz lasing is the accumulation of a 
large population of THz photons (n c ^> 1), for which a strong 
nonlinearity and a good THz cavity are required. Further- 
more, this should also lead to a coherence buildup (g K c ' — )> 1). 
Assuming that other destructive issues (as pure dephasing) 
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are negligible in the system, the efficiency of such coher- 
ence buildup is further enhanced by the irreversibility of the 
polariton-to-photon conversion ifTTIl . The conversion should 
happen effectively mainly in the polariton-to-photon sense so 
that the formation of correlations is not disrupted by the in- 
verse process. However, if, after one upper polariton has 
transformed into a THz photon plus a lower polariton, these 
two last remain for a long time in the cavity, the inverse pro- 
cess is more likely to happen. Given that, at the same time, 
we wish the THz population to remain high, the only way to 
suppress the photon-to-polariton process is that lower polari- 
tons do not remain for long times in the cavity, that is, we 
need Y^ Yc- Therefore, we have chosen G = y and y c = 0. ly, 
which allows for an accumulation of photons faster than of po- 
laritons, as shown in Fig.[2ja). For the moment, we have also 
neglected the competing process through the exciton reservoir 
(Yab — 0) mat does not produce THz radiation and clearly hin- 
ders lasing. In these conditions, Fig. [5Ja) shows a nonlinear 
increase of the THz population (circles) which is accompanied 

by the expected g^ —> 1 in Fig. [^[b). Note that there is not 
a sharp threshold for the lasing emission. The THz photons 
gradually and monotonously increase and acquire coherence. 
One can in any case define the threshold as the pump where 
n c = 1, that is, P w 0.22y in this case. 

To link the present full-quantum master equation approach 
with the kinetic model (8), we derive from Eq. ([T]) the steady 
state equations for the average populations, that are of the type 
of a Boltzmann source and sink rate equation, i.e., 

n x = Pf/lf (3) 



for x = a,b,c. This provides us with the effective pumping 
and intensity decay rates for the three modes: 
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The semiclassical rate equations of Ref. are recovered by 
setting all cross correlation functions to unity. Solving exactly 
these rate equations, one obtains approximated formulas for 
the populations that, in turn, provide analytical expressions 
for all the effective parameters above. The expressions are too 
lengthy to be given here but we still note them as ^ ate , P™ iQ 
and FJ. ate for later comparison. 

In our full quantum derivation, cross correlation functions 
are computed numerically and self-consistently and can be, in 
general, larger (resp. smaller) than unity, showing bunching 
(resp. antibunching) of joint emission in the modes x and y. 
This is shown in Fig. [2jc), where cross emission between 
upper polaritons and THz photons is antibunched (squares, 

(2) 

gac < 1), as the destruction of the first one implies the cre- 
ation of the second ones. This also tends to be the case be- 
tween upper and lower polaritons, however, in this case due to 
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FIG. 3. Relevant THz quantities in the steady state as a function of 

(2) 

pump: n c (first line), g K c (second line) and T c (third line). The first 
column (plots a, b, c) corresponds to fixing y c = O.ly and increasing 
Yah in the sense indicated by the arrow (from to 0.5 in steps of 0.05). 
The second column (plots d, e, f) corresponds to fixing Yab = an d 
increasing y c in the sense indicated by the arrow (from 0.1 to 0.9 in 
steps of 0.1). 

the large difference <C ft c , they have become independent 

(2) 

(rhombus, g y ab ^ 1). The cross emission of lower polaritons 

(2) 

and THz photons is bunched (circles, g K bc J > 1) for the same 
reason: they are produced at the same time and their emission 
tends to be simultaneous. At very high pumps all the emis- 

(2) 

sions become statistically independent (gx y J — >> 1). 

Another typical manifestation of lasing is line narrowing 
of the luminescence emission, that is, the increasing lifetime 
of the photons emitted by the lasing mode: in the incoher- 
ent regime, the linewidth corresponds to the inverse lifetime 
of a single (incoherent) photon, whereas in the lasing regime, 
THz photons form a collective state with a lifetime affected by 
the total population. The effective inverse lifetime of the THz 
mode is given approximately by Eq. ( |4f| ) in the limit of inde- 
pendent classical modes. T^ ff decreases appreciably when G 

(2) (2) 

is large and gL flfe ~ gac n a decreases with pumping. In the 

case of Fig. [2J where g a 2 ) « 1 and n a , are small, the main 
element that determines the linewidth is the cross correlation 

(2) 

function g K b J , which is large (bunched) but strongly decreasing 
with pump. 

To check that this process is accompanied by the expected 
line narrowing of the emission, predicted by Eq. ( |4f| ) in the 
regime of Fig. [2j we compute the photoluminescence spec- 
trum of THz emission. It is defined, in the steady state, as 
S(co) oc 9i J °°(c t (0)c(T))^ a)T (iT and is computed (with the 
quantum regression formula) from the two-time correlator: 

(c+(0)c(T)>= £ V~rW [n , mA {x) (5) 
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in terms of W^ m ^ , that follow the same equation as the off- 
diagonal density matrix elements (n,m,r\p \n,m,r—l), with 
initial values W^ m ^ (0) = V^^[^,m,r] • The required equations 
read: 

= - { 7(n + m) + Yc (r - i ) + + 1 ) } 
+ y{(n+ 1)^+1^] + (m+ l)*[„^+i,r]} 

+ 7ofe { (l + l)^^+l,m-l,r-l] - n ( m + 1)>^[ W ] } 
+ G(n + l)m| VK r - 1 )^[n+l,w-l,r-l] - ( r ~ \W[n,m,r] } 

+ Gn(m+l)| v / r(r+l)^_ 1?m+1?r+1 ]-(r+^)^ m?r ]|. 

(6) 

The THz spectrum of emission consists of only one peak. The 
associated linewidth is extracted and plotted in Fig.[2jd). We 
find that, although not quantitatively exact, the trend (narrow- 
ing) follows remarkably that predicted by the effective THz 
intensity decay rate, Eq. pf) . 

Finally, in Fig.|3j we analyse how two key factors influence 
the quality of the THz lasing. Still as a function of pump, we 

(2) 

show how THz population (first row), coherence g K c 1 (second 
row) and its spectral linewidth, T c , (third row) are spoiled as, 
on the one hand, incoherent relaxation from upper to lower 
polaritons through the reservoirs (i.e., not through the THz 
mode) is increased (left column), and, on the other hand, as 
the THz cavity Q factor is decreased (right column). The 
sharpest feature to tell when lasing disappears is the slope of 
the linewidth curve, that becomes positive at all pumps. This 
happens at Yab ~ 0.2 y in the first case and at y c « 0.7y in the 
second. These are the minimum values required to achieve 
the THz lasing with all its features. The analytical approxi- 
mated expression for the linewidth, T^ ate , is only valid for a 
bad THz cavity, when the system is in the "weak coupling 
regime" (y c > G, y) with the THz mode (see Fig. [4]). Here, 
the rate equations, that neglect quantum correlations, give a 
good description of the dynamics. With better systems able 
to enter the lasing regime, one must solve the full quantum 
problem, including one- and two-time correlators, in order to 
retain important qualitative features of the linewidth such as 
line narrowing. 



To conclude, our analysis shows that although an efficient 
THz lasing is indeed possible in optically pumped polariton 
systems, it requires a special care in the design of the struc- 
ture. One must ensure a good THz cavity and small deple- 
tion of the k = state of the upper polariton mode due to the 
phonon-assisted relaxation of polaritons from this state to the 
states of the lower polariton mode characterised by large in- 
plane wavevectors k. The latter obstacle is the most serious as 
upper polariton depletion is usually very fast in planar cavi- 
ties. This can be circumvented by the use of pillar microcavi- 
ties ifTUl , where a full three-dimensional photonic confinement 
is achieved. In pillars, the spectrum of exciton-polaritons con- 
sists of discrete states, so that there is no reservoir of states 
0.5 R 
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FIG. 4. Difference between the numerically computed THz emis- 
sion linewidth (Fig. |3j:,f) and the analytical approximated formula 
obtained from the rate equations: 100(r£ ate — r c )/r c . Pumping is 
fixed at P = 0.451 y. 



with large wavevectors. Optical transitions between exciton 
polariton states in the pillars are usually forbidden, but if one 
of these states is brought into resonance with the elhh2 or 
e2hhl exciton state, the transition between this state and po- 
lariton states produced by coupling of elhhl excitons with 
light becomes possible with emission of a THz photon. Also, 
a pillar microcavity is certainly a more appropriate structure 
to be fitted in still another cavity, the THz one, needed for THz 
lasing. 
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